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We discuss theoretically the behavior of the velocity au- 
tocorrelation function in the Dissipative Particle Dynamics 
model. Two dynamical regimes are identified depending on 
the dimensionless model parameters. For low frictions a mean 
field behavior is observed in which the kinetic theory for DPD 
provides good predictions. For high frictions, collective hydro- 
dynamic effects are dominant. We have performed numerical 
simulations that validate the theory presented. 



I. INTRODUCTION 

The Dissipative Particle Dynamics model allows to 
simulate hydrodynamics at mesoscopic scales in which 
thermal fluctuations are important jf]||]. For this rea- 
son, it appears as a good simulation technique for the 
study of complex fluids like polymer or colloidal suspen- 
sions where both hydrodynamic interactions and Brown- 
ian motion are important . Being an off lattice tech- 
nique, it does not suffer from the restrictions imposed by 
the lattice as in Lattice Gas Cellular Automata or the 
Lattice Boltzmann approach and it is much more flexible 
for modeling. 

Even though the technique has a well sounded theoreti- 
cal support and many applications have been undertaken, 
there is at present no systematic study of the region of 
parameters much suitable for simulation of particular hy- 
drodynamic problems. In addition, recent simulations 
H have shown deviations from the transport coefficients 
predicted by the kinetic theory developed by Marsh et al. 
[]9| . The two approximations involved in this kinetic the- 
ory are the small gradient expansion around local equilib- 
rium and the molecular chaos hypothesis. However, it is 
difficult to investigate the origin of the discrepancies be- 
tween theory and simulations within the kinetic theory 
context. The theory just produces the explicit expres- 
sions for the transport coefficients with no hint about its 
range of validity. It has been suggested that it is pre- 
cisely in the region of parameters where kinetic theory 
fails where it is more sensible to conduct simulations that 
reproduce hydrodynamic behavior J8|. 

We shed some light into the problem by presenting a 
theory that allows to compute the velocity autocorrela- 
tion function (vaf) of the dissipative particles. The the- 
ory is based on the physical picture in which the DPD 
particles are regarded as Brownian-like particles moving 
in an environment created by the rest of DPD particles. 
Strictly speaking, though, they are not Brownian parti- 



cles because the total moment of the system is conserved. 
This approach was introduced by Groot and Warren as a 
way of computing the diffusion coefficient |j] . We identify 
the basic dimensionless parameters which allow to clas- 
sify and discuss the dynamical regimes displayed by the 
model. By assuming that the environment of DPD parti- 
cle behaves hydrodynamically, it is possible to obtain an 
explicit analytical expression for the velocity autocorrela- 
tion function. Finally, we present numerical simulations 
that allow to validate the hypothesis made in the theory. 



II. THE DPD MODEL 

The stochastic differential equations (SDE) that gov- 
ern the position and velocity of the «-th particle of 
mass m in DPD are given by 
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where the following quantities are denned 
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In order to compare with the kinetic theory in Ref. ||, 
it is assumed that the usual conservative force is not 
present. The noise amplitude a is given by the detailed 
balance condition 



a = (2kBTjm) 



1/2 



(3) 



where T is the temperature of the equilibrium state to- 
wards which the system relax (if the boundary conditions 
allow for it) and ks is Boltzmann's constant. Finally, 
dWij = dWji are independent increments of the Wiener 
process that obey the ltd calculus rule 

dWijdWi>j> = (5u'5jj> + Sij'Sji')dt, 



(4) 



i.e., dWij is an infinitesimal of order 1/2 ]l0| . The dimen- 
sionless weight function uj{r) is normalized according to 
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where n is the number density of the system. In this 
paper we will work in two spatial dimensions (2D) and 
select the following weight function with range r c , 
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if r < r c and zero if r > r c . 

We discuss now which are the fundamental parame- 
ters for the DPD model. By appropriate choice of units 
of mass, time and space it is always possible to reduce 
the number of relevant parameters of the model. It is 
obvious that the dynamical regimes are independent of 
the units used, and will depend on dimensionless pa- 
rameters only. There are six parameters in the model: 
m, 7, r c , ksT 1 A, L, where A is the average distance be- 
tween particles, related to the number density n of parti- 
cles by A = n-^ d , dis the space dimension, and L is the 
box size (or any other boundary length scale). From these 
six parameters we can form three dimensionless parame- 
ters. By defining the thermal velocity vt = (hsT/m) 1 / 2 
we select 



di>T dr-y 
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The physical meaning of these parameters is as follows: 
tt is the time taken by a particle moving at the thermal 
velocity to move a distance r Cl whereas r 7 = 7" 1 is the 
time associated to the friction. Therefore, the dimension- 
less friction fl is the ratio of these two time scales. On the 
other hand, s is the overlapping between particles which 
is related to the number of particles that arc within the 
range of interaction (the action sphere) of a given one. 
Finally, /i is the dimensionless box length. These dimen- 
sionless parameters i7,s,// fix the dynamical regimes of 
the model. 



III. VELOCITY AUTOCORRELATION 
FUNCTION 

The velocity equation in ([In can be written in the form 



^w(r y )( 
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where the random force is Fidt = ^ ( r ij) e ijdWij . 
We have introduced in ^ the environment velocity by 

Vj (t) = d ^ u( r v ) e ij e ij ' v i • ( 9 ) 



This velocity is a weighted average of the velocities of 
the neighboring particles of particle i. Next, we observe 
that the factor of in the right hand side of (||) can be 
written as 

v ^tj(r ii )eyey = / dru^r, -r)-^ — ^ — -rn(r,t), 
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where we have introduced the microscopic density field 
n(r,i) = J2j^i 3( r j ~ r )- ^ we assume that the density 
field is constant with value n (which will be confirmed by 
the results obtained later) then we may approximate 



druj(ri — r) 
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The last equality is obtained by noting that the integral 
is an isotropic second order tensor, which must be pro- 
portional to the identity (the constant of proportionality 
can be obtained by taking the trace of the integral and 
using the normalization (]q)). 

After using Eqns. (|Io| ) and ([ll]) in Eqn. @ one obtains 



dr., =v t dt, 

dv i = -'^[v i -V i ]dt+—dt. 
a m 



(12) 



We observe that the DPD particles behave similarly to 
Brownian particles but in a systematic velocity field de- 
termined by the rest of its neigbouring particles. The 
stochastic properties of the random force are not exactly 
those of a Brownian particle because the total momen- 
tum of the system is conserved but for the rest of the 
development they are irrelevant. 
The formal solution of Eqn. (|l^) is 

Vi(t) =cxp{- 7 i/d}v l (0) 



dt'exp{-j(t-t')/d} 
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By multiplying this equation by Vi(0) and Vi(0) and av- 
eraging one obtains a set of equations for the velocity 
autocorrelation function at equilibrium 
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+ l f dt' exp{- 7 (t - OM^ViCtO-VitO)), 
^<v i (t).V i (0)) = 

\ J* dt'exp{- 7 (* - )-V,(0)>, 
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where use has been made of the fact that the random 
force is not correlated with the velocity at present and 
past times and the property (Vj(0) -Vj(0)) = (which 
can be checked from the definition (||)). Substitution 



of the second equation in (14) into the first one leads 
to an expression that relates the particle vaf with the 
environment vaf, this is 

i< Vi (i).v,(0))=«p{-7t/d} — 
a m 

+ {if / dt ' {t ~ t>] exp {~^ - OM^v^tO-v^o)). 

(15) 

The second term in the right hand side represents collec- 
tive effects. When this term is negligible we say that a 
mean field approximation is valid, in which the velocity 
autocorrelation function decays exponentially. The rea- 
son for the name "mean field" comes form the observation 
that Eqn. (12), in which the average value (Vj) = is 
used instead of the instantaneous value Vj, produces an 
exponential decay of the velocity autocorrelation func- 
tion. In the Appendix of Ref. [|| it was computed the 
velocity autocorrelation function and the diffusion co- 
efficient of the DPD particles by using this mean field 
approximation. 

The solution (|l5|) is still formal because we do not know 
explicitly the form of the correlation of the environment 
velocity (which will be given in the next section) . Never- 
theless, it is possible to extract useful information from 
this expression. This is most conveniently done by taking 
dimensionless variables. Let t be the dimensionless time 
tvT /r c , this is, the time expressed in units in which r c — 1 
and vt = 1, and v = v/vt a dimensionless velocity. In 
these units, Eqns. (H4) take the form 



1 



yi(t)-vi(0)) = exp{-fii} 



In order to complete the picture we need to evaluate the 
overall magnitude of the correlation of the environment 
velocity. This is determined by the value at the origin 
which is computed in the Appendix I with the result 



+ il / dt exp{-n(?-f )}-(Vi(T)-Vi(Q)), 

Jo d 
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and Eqn. ( |l5| ) takes the form 

-<Vi(?)-Vi(0)) =exp{-Ol} 

ft 1 _ _ 

+n 2 / dt'(t-t')exp{-n(t-t')}-(Vi(t')-V t (0)) (17) 
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For later notational convenience we introduce 



c(t) = -<Vj(i)-Vj(0)>, 



c(i) = I< V j(t).Vj(o)>. 



i<V i( 0).V,(0)> = ^. 



(19) 



We observe that the magnitude of this correlation de- 
creases with the overlapping coefficient. This is phys- 
ically meaningful because the environment velocity is a 
weighted average of the velocities of the particles that are 
within an action sphere. These velocities are distributed 
at random and, therefore, if there are many particles 
within an action sphere, the average will be proportion- 
ally smaller. 

Now, several qualitative predictions concerning the dif- 
ferent dynamical regimes can be made from expression 
(|l6| ) or (17). For fixed 0, the large overlapping s limit 
produces a small contribution from the collective part 
and the velocity correlation function decays in an expo- 
nential way. For fixed overlapping s, when f2 is small 
(in the limit of small friction or high temperature) the 
behavior of the vaf is again exponential. In the opposite 
regime of large il, the exponential contribution decays in 
a very short time and the main contribution for times 
larger than is given by the collective term. Actually, 
in the limit f2 — > oo the exponential memory function 
acts as a delta function and for times larger than f2 _1 
one obtains 

(vj(t)-Vj(O)) « (Vj(I)-Vj(O)) « (Vi®-Vi(0)). (20) 



The physical meaning of these expressions (|20|) is also 
clear. When the friction is high, in a very short time 
the velocity of a given particle is slaved by the average 
velocity of its environment. 



IV. HYDRODYNAMIC HYPOTHESIS 

The environment velocity Vj defined in Eqn. (^|) can 
be rewritten as an average over the action sphere of the 
microscopic velocity field, this is, 



(16) 



Vi(t) = d / drwfo-r)- 



r n(r,i)v(r,t), 



where 



n(r,*)v(r,i) = ']Tv j 5(r :j - r). 



(21) 



(22) 



(18) 



The velocity field v(r, t) obeys the equations of hydrody- 
namics when its characteristic length scale is much larger 
than the interparticle distance. We expect that the av- 
erage involved in Vj will be dominated by the hydro- 
dynamic modes whenever the range of interaction r c is 
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much larger than the interparticle distance A (i.e., large 
overlapping s). In this section we present a calculation 
of the correlation function of the environment velocity 
which is based on this hydrodynamic argument. 

From Eqn. (pl|), the environment velocity is conve- 
niently expressed in terms of the Fourier components of 
the velocity field v(k, t), this is 

/dk 
— ^(k^v^exp^kri}, (23) 

where we have introduced the second order tensor 



;(k) = J dru>(r)rr cxp{— ikr}. 



(24) 



The explicit form of this tensor when the weight function 
is given by Eqn. (^|) is given in the Appendix II. The 
environment velocity correlation function is given by 



(V 4 (0)-V. t (t)) = (dn 



dk dk' 



w(k)w(k') 



(2tt) 2 (2tt) 2 
x (exp{ l (kr,(0) +k' ri (t))}v(k,0)v(k',t)>. 

(25) 

We will assume that the position of particle i is weakly 
correlated with the velocity field v(k, t) in such a way 
that we can approximate 



(exp{i(kr i (0) + k' ri (t))}v(k, 0)v(k', *)) 

« (exp{i(kr i (0)+k / r i (t))}) (v(k, 0)v(k', t)). 



(26) 



We further assume that the correlation of the velocity 
field is given by the linear hydrodynamics result 



<v(k',0)v T (M)> 



(27T) 2 5(k + k') 



x[exp{-uk 2 t}(l - kk) +exp{-rfc 2 t}cosfccikk]. (27) 

Here, v is the kinematic viscosity, T is the sound absorp- 
tion coefficient, and c is the sound speed of the DPD fluid. 
The correlation function ( |27| ) is different from zero only 
when k = — k'. The first average in the right hand side 
of ( p(| ) is, therefore, given by the incoherent intermedi- 
ate scattering function F s (k, t) = (exp{£k(ri(0) — 
|TT|| . By further assuming a hydrodynamic behavior for 
this function we obtain [fill 



F s (k,t) = ex_p{-Dk 2 t}, 



(28) 



where D is the self-diffusion coefficient of the DPD par- 
ticles. 

The final hydrodynamic expression for the environ- 
ment velocitycorrelation function is found by substitu- 
tion of ©,@,(fl, into © 
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where the following functions are defined 
J dka 2 (k) exp{— xk 2 } 



$(ar) 



Jdka 2 {k) 
J dk(a(k) + b(k)) 2 exp{— yk 2 } cosfcz 
Jdk(a(k) + b{k)) 2 ' 



(30) 



which satisfy $(0) = 1,*(0, 0) = 1, with the definitions 
for a(k), b(k) given in Appendix II. 

In this way, the environment velocity correlation func- 
tion is explicitly given in terms of hydrodynamic fluid 
properties, i.e., the transport coefficients of the DPD 
fluid. The prediction is not complete until particular 
values for these transport coefficients are provided. One 
possibility is to measure these transport coefficients in a 
simulation. Another possibility, which is the one we fol- 
low here, is to use the values for v, D, T, c provided by ki- 
netic theory JsJ . This has the advantage that one knows 
the explicit dependence of the transport coefficients in 
terms of the dimensionless parameters f2, s. Even though 
for some dynamical regimes the kinetic theory predic- 
tions are not in exact quantitative agreement with the 
measured transport coefficients ||,[l2j the discrepancies 
between theoretical and simulation transport coefficients 
are small. 

The kinetic theory results are H, ]13|| 
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The sound absorption coefficient is given by T = 2v+^V[, 
for the case that the equation of state is that of the ideal 
gasQ. 

We are now in position to write the hydrodynamic pre- 
diction (pS[) i n terms of dimensionless variables. Substi- 
tution of (|3l|) with the expression (^) into (^) leads to 
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(32) 



We observe that the time scale of the environment ve- 
locity correlation function is determined by fi, while its 
amplitude is determined by s. In Fig. |Vj we show the 
theoretical prediction ([32]) for a particular value of s and 
three different values of O. An algebraic dependence as 
t _1 is observed for very long times which is the celebrated 
long time tail t~ d / 2 arising from the diffusive character 
of the shear mode. 



4 



io- 



0.01 



0.1 



1 



10 



100 1000 



t 



FIG. 1. Theoretical hydro dynamic prediction for the en- 
vironment velocity correlation function C(t) for a value of 
s = 2.82 and three different values for Q. Lowest curve is for 

= 25, middle curve is for Q = 8.3, and upper curve is for 
fl = 0.5. The algebraic long time tail behavior appears for 
very long times, typically when the correlation has decayed 
three orders of magnitude from its initial value. Dotted line 
is a curve proportional to t~ x to guide the eye. 



V. SIMULATION RESULTS 

We have simulated Eqns. ([!]) in two spatial dimensions 
in a system with periodic boundary conditions. The ve- 
locity autocorrelation function of the DPD particles and 
also the environment velocity autocorrelation function 
have been computed at equilibrium. 

In order to derive Eqn. ( p"7| ) we made the approxi- 
mation that the density field is approximately constant. 
We check now that this assumption was reasonable by 
computing in a simulation the environment velocity cor- 
relation C(i), evaluating numerically the integral term in 
Eqn. (|17|) , and adding up the first exponential term in 
Eqn. fllTf). The result is the dotted line in Fig. ||. Also 
shown is the result for the velocity correlation function 
c(t) obtained directly from the simulation (solid line). We 
see that both results are in quite good agreement giving 
confidence on Eqn. ( |l7| ) as a sounded starting point for 
theoretical analysis. Similar good agreement is obtained 
for all the values of £1 studied (17 = 0.5, 8.3, 25) 
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FIG. 2. Comparison between the simulation result for c(t) 
(solid line) and the value given by Eqn. (j^) (dotted line) us- 
ing the simulation result for the environment velocity corre- 
lation function C(t). Parameters: s — 2.82, Q — 25, /i = 10.0 

For a given value of 51 = 25 we plot in Fig. [y] the 
value of s 2 C(t), for different values of s. According to the 
hydrodynamic prediction ( |32| ) this curve should be inde- 
pendent of the overlapping coefficient s (dotted curve in 
Fig. [v]) . We observe that the simulation results converge 
towards the theoretical prediction only when the over- 
lapping coefficient is sufficiently large. This is expected 
from the fact that hydrodynamic behavior appears only 
on length scales that involve a relatively large number of 
particles. 



s 2 C{t) 




FIG. 3. Simulation results for the velocity correlation func- 
tion s 2 C(t) for a fixed value of Q, — 25 and different values 
of the overlapping coefficient s = 1.5,2.18,3.2 (upper curve 
correspond to lower s). Dotted line is the hydrodynamic pre- 
diction (^). As the overlapping increases, the simulation 
results converge towards the theoretical prediction. 



We have investigated also the effect of the finite system 
size. We plot in Fig. ^ the environment velocity correla- 
tion function for different system sizes, while keeping the 
rest of parameters constant (s = 2.82, = 25). We ob- 
serve a large discrepancy between the hydrodynamic pre- 
diction and the simulation results when the box is small. 



5 



This discrepancy appears at large times and is reduced 
when the system size increases. This effect can be un- 
derstood as an self-interaction through periodic bound- 
ary conditions. According to Onsager hypothesis for the 
regression of fluctuations, the correlation of the veloc- 
ity fluctuations decays in essentially the same way as a 
localized macroscopic velocity perturbation. A pertur- 
bation in the velocity field in the form of an impulse at 
the origin decays through the formation of a vortex ring 
whose center evolves diffusively as y/i [Q. In a system 
with periodic boundary conditions the vortex ring is rein- 
troduced into the system producing an effective slowing 
down of the velocity field at the origin. Obviously, the 
larger is the system size, the latter the effect appears and 
also the smaller is the amplitude of the effect. 
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FIG. 4. Simulation results for the velocity correlation func- 
tion C(t) for fixed value of s = 2.82, fi — 25 and different 
values of the box size [i = 7.9,11.2,15.8,35.5 (lower curve 
correspond to lower fi). Thicker line is the hydrodynamic 
prediction (|32|). 



Next, we study the effect of fi on the velocity autocor- 
relation function of the DPD particles. In Fig. [Vjwe show 
the velocity correlation function for given s = 2.82, [i = 
1 0.0 and three different values of fl = 0.5, 8.3, 25. We also 
plot the corresponding exponential terms in Eqn. (17). 



For small Q (small friction or high temperature) the de- 
cay of the vaf is very accurately given by the exponential 
term. As long as fi increases, discrepancies from the ex- 
ponential behavior are observed. This discrepancies are 
due to the effect of the collective term in Eqn. (|l7|). 
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FIG. 5. Simulation results for the velocity correlation func- 
tion c(I) for fixed value of s = 2.82, fi = 10.0 and different 
values of fi = 0.5 (circles), fi = 8.3 (triangles) and fi = 25 
(diamonds). Solid lines are the term exp{— fit} in Eqn. (|l7|). 

In Fig. it is shown that for large CI, at large times the 
environment and particle velocities coincide according to 
Eqn. ©. 




FIG. 6. Simulation results for the velocity correlation func- 
tion c(t) (thin line) and the environment velocity correlation 
function CCt) (thick line). Also shown is the term exp{— fit} 
in Eqn. rtTt). Here, s = 2.82, fj, = 10.0, fi = 25. 



VI. SUMMARY AND DISCUSSION 

We have presented a theory for Dissipative Particle 
Dynamics that allows to understand the different dy- 
namical regimes displayed by the model. The theory is 
based on the physical picture that DPD particles behave 
like Brownian particles in a non-equilibrium environment 
due to the rest of DPD particles. An explicit expres- 
sion for the velocity autocorrelation function is derived in 
which the Brownian exponential behavior is corrected by 
the presence of collective effects. By using dimensionless 
variables it is possible to asses the range of parameters 
for which the collective effects are important. Three di- 
mensionless parameters appear in the model, s, CI, fi, and 



G 



they characterize the different dynamical regimes in the 
system. The relevance of precisely these dimensionless 
groups is motivated by the theory presented (in Ref. ||] 
other dimensionless groups are introduced). 

Two dynamical regimes are identified, the mean field 
regime and the collective behavior regime. The transition 
between both is governed essentially by the dimensionless 
friction CI with important effects of s and /z. Mean field 
behavior appears for small friction f2 or large overlapping 
s. This is physically reasonable: For small friction, the 
dynamics of the environment of a given particle hardly 
affects the behavior of this particle. For large overlap- 
ping the non-equilibrium collective effects are smeared 
out over large regions. The mean field approximation is 
closely related to the molecular chaos assumption made 
in kinetic theory. Actually, it is possible to compute 
the diffusion coefficient of the DPD particles by using 
the mean field approximation for the vaf into the usual 
Green-Kubo formula for the diffusion coefficient || . The 
result is precisely the prediction for the diffusion coeffi- 
cient given by kinetic theory [3 (with due account of the 
normalization given in Eqn. (E)), i.e. 



D = 



1 , , . ... d knT 
5 <vi(t)-v i (0)) = -^-. 



(33) 



When the dimensionless friction Q is high, the vaf does 
not decay in an exponential way because it is dominated 
by the collective dynamics. We have presented a theo- 
retical prediction for the collective part of the vaf by as- 
suming that the correlation function of the environment 
velocity reflects an underlying hydrodynamic behavior. 
Such a behavior is expected (and observed in the simula- 
tions) when the overlapping coefficient s is large enough. 
In this case, the collective effect is small but well de- 
scribed by hydrodynamics. 

The fact that hydrodynamics governs the dynamics 
of the velocity of the particles has two important con- 
sequences. The first one is that finite size effects are 
important. Hydrodynamic self- interaction through the 
periodic boundary conditions exists and large box sizes 
must be considered in order to render this effect small. 
This self-interaction occurs not only through the sound 
mode but also through the shear mode jl5| . The second 
consequence is the appearance of the celebrated long- 
time tails in the vaf. These algebraic tails occur at very 
long times (for which the vaf has decayed to a factor 
10 -3 from its original value) and are difficult to observe 
in our simulations due to the statistical noise. Never- 
theless, we have provided sufficient numerical evidence 
for the hydrodynamic behavior at smaller times and we 
expect the presence of long time tails at large times for 
sufficiently large system sizes. The i" 1 dependence of the 
vaf, when introduced into the Green-Kubo formula j3^ ) 
leads to a logarithmically divergent diffusion coefficient. 
It is apparent that this "small" divergence will be hardly 
observable in any simulation with a finite box size. 

In the regimes for which collective effects are impor- 
tant, we expect that deviations from the predictions of 



kinetic theory occur not only for the diffusion coefficient, 
but also for the rest of the transport coefficients of the 
DPD fluid ||. It is actually possible to use the theory 
presented in this paper in order to compute the trans- 
port coefficients for the DPD fluid expressed in the form 
of Green-Kubo formulae jlqj . One obtains then a set 
of recursive equations in which the transport coefficients 
are expressed in terms of the transport coefficients used 
in the hydrodynamic assumption. We do not follow this 
rather cumbersome numerical procedure here because it 
does not provide any new physical insight. 

Dissipative Particle Dynamics is designed to simulate 
hydrodynamic problems. Actually, one would like to have 
the DPD particles moving in such a way that they fol- 
low accurately the flow field intended to be modeled. We 
see that this will happen whenever the friction is suffi- 
ciently large (in this case the velocity of a DPD particle 
is slaved by its environment) and the overlapping is suffi- 
ciently large (in such a way that the environment velocity 
moves hydrodynamically) . In this regime, the dynamics 
of the particles is mainly collective and kinetic theory 
gives inaccurate values for the transport coefficients . 
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APPENDIX I 

We derive the equilibrium value for the environment 
velocity correlation function given in Eqn. From 
the definition of V; in Eqn. (^|) we calculate the mean 
value by using the canonical ensemble 

^(V;(0)-Vi(0)) = d((^2 UijGijeij-Vj)(^2 uj lk e ik e ik -v k )) 



dk B T(N- 1) j dT t dT 3 uj{ 
m J dYidYj 

dk B TN~l 3 

m N lixr^n 
dk B T 3 

m 2its 2 ' 



(34) 



for large N (number of DPD particles in the system) and 
n the number density. 



APPENDIX II 

On simple symmetry grounds, the tensor u>(k) has the 
form 
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ui(k) = J druj(r)rr cxp{— ikr} = a(kr c )l + 6(fcr c )kk 

(35) 

In order to calculate the functions a(kr c ) and b{kr c ) we 
double contract w(k) with the dyadic kk and also take 
its trace. This leads to 



k T cj(fc)k = a(kr c ) + b(kr c ) 
= 2tt 



druj(r) ( — - — rJiikr 



tr(w(fc)) = 2a(kr c ) + b(kr c ) 

poo 

= 2tt ruo{r)J {kr). 
Jo 



(36) 



The integrals are given in terms of generalized hyperge- 
ometric functions p F q {&, b, z} |l7]] and Bessel functions 



p F q {a, b, z} = ^2 



(ai)fc...(a p )fe z k 



k ^ {h)k--(b q )k k\ 
where (m)k = m(rn + l)...(m + k — 1). More precisely 



(37) 



a(Ar c ) + b(kr c ) 



■i.*«§).ftf),-^> 



G 



(> 



:Jo{kr c 
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{kr c f 



■Ji{kr c ) 



2a(kr c ) + b(kr c ) 
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SMkr c ) 2 3 5 (fcr c ) 2 
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The solution of this system of two equations provides the 
values for a(kr c ) and b{kr c ). 
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